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Abstract 

The paper describes the derivation of the Kohn-Sham equations for a nanowire with direct 
current. A value of the electron current enters the problem as an input via a subsidiary condition 
imposed by pointwise Lagrange multiplier. Using the constrained minimization of the Hohenberg- 
Kohn energy functional, we derive a set of self-consistent equations for current carrying orbitals of 
the molecular wire. 
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I. INTRODUCTION 

The goal of building real electronic devices from individual molecules Q has spurred 
experimental and theoretical studies of molecular conductivities. Measurements of current- 
voltage (I-V) characteristics have been reported for single molecules but mechanisms 
for electron transmissions through molecular interfaces remain largely unexplored and un- 
explained 0, 1^ . It appears that the miniaturization of microelectronic devices shrinks the 
size of the system to a level when the simulations of transport characteristics of molec- 
ular electronic devices can be performed by ab initio electronic structure methods (with 
single benzene- 1,4-dithiolate molecule as the molecular wire, the contact wire-contact "su- 
permolecule" consists of just 32 atoms). However, the main principle hurdle for the appli- 
cation of the first principle electronic structure methods is not the size of the system but 
conceptually unsolved problems of modeling of the electronic current and interpretation of 
the applied voltage for correlated electronic systems. In traditional methods, which are 
based on the Landauer theory Q], the molecule and the edge part of the metal contacts can 
be in principle treated ab initio, e.g. in the density functional theory (DFT). But then to 
induce the electronic current it is assumed that the system is attached to two equilibrium 



electronic reservoirs with the different chemical potentials j7(. Being computationally not 
very efficient, since it requires to perform the calculations over the whole coordinate space, 
these traditional approaches are also pivotally relied on the electronic reservoir concept and 
often use ad hoc reservoir-molecule interaction as an essential input parameter. 

We have previously shown that the reservoir model can be conceptually circumvented in 
electronic transport calculations if the electronic current is used as an input information in- 
stead the applied voltage bias 0, 1^ . In our approach the electronic current enters a quantum 
variational problem via the corresponding Lagrange multiplier. In this paper, we continue 
development of our method and discuss detailed derivation of the Kohn-Sham equations for 
a nanowire with steady current. We begin by describing the variational principle for an 
inhomogeneous electron gas with steady current. We then derive the Kohn-Sham equations 
with steady current. 



2 



II. VARIATIONAL PRINCIPLE 



We start our derivation by defining density and current density of an infiomogeneous, 
interacting A^-electron gas in an external scalar potential v{r) as (a.u. are used tlirougfiout 
tlie paper): 

N 

p(r) = iV / vl/(r,r2,...;r^)v[/*(r,r2,...,r^)n^r, , (1) 

N 

j(r) = N Im{^*(r,r2,...,r^)V^(r,r2,...,r^)}n^r, , (2) 

i=2 

wliere \1/ is the many-electron wave-function. For simplicity of notation we do not write 
explicitly the spin indeces but our derivation can be readily extended to spin-polarized 
systems. 

Next, we introduce a subsidiary condition for the current density distribution by speci- 



fying the quantum wire constraint on the current 



dydzj^{r) = I{x), (3) 

where I{x) is a given input value of the steady current through the electronic device. With 
this constraint the current is aligned along the x-axis and the net current flow across a cross 
section J dydzjx{r) is constrained, and this quantity is available experimentally. 

For an inhomogeneous, interacting electron gas in scalar external potential v{r), the 



ground state energy is given by the Hohenberg-Kohn energy functional 
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EhkIp] = TM + I dTv{r)p{r) + ^ / / ^rdr'^^^^ + E^M ■ (4) 

The first term, Tq, is the kinetic energy functional of noninteracting electrons with given 
density and current density distribution. The second and third terms, are the interaction 
energy with external potential and the electrostatic interactions. The last term, E^j-dp], 
is the exchange and correlation energy functional. The current dependent correction to 
the standard exchange-correlation energy functional is O(j^) A dependence of the 

exchange-correlation energy functional upon j(r) can be readily included in our approach, 
yielding the additional exchange-correlation vector potential 



A^iP,m = (5) 



which modifies the kinetic energy operator in the final Kohn-Sham equation in the same 
way as an external vector potential does. The current dependent correction to the standard 
exchange-correlation energy functional is 0{j^) Molecular and atomic wires can sustain 
a current of about 10~^a.u. for typical experiments on molecular conduction 2]. Therefore 
the direct current-dependent corrections to the total electronic energy are only of the order 
of 10~^a.u. and can be ignored in most molecular transport calculations. 

Likewise our derivation of the Schrodinger equation for current carrying states jol, we 
use the constrained variational principle in our derivations. The problem is to minimize 
Hohenberg-Kohn energy functional subject to the constraint on the current density Q. 
The constraint is explicitly achieved by introduction of an auxiliary functional 



n[^] = Ehk[p] - j dxa{x) (^J dydzj^iv) - /(x)^ . 



(6) 



The first term giving the total energy is standard in the variational formulation of the DFT. 
The second term with the Lagrange multipliers a{x) has been introduced to impose the 
subsidiary constraint for the current density 



III. KOHN-SHAM EQUATIONS WITH DIRECT CURRENT 



Following the Kohn and Sham approach we introduce a reference fermion system 



with orthonormal single-particle orbitals i/jiir) and occupation numbers rii to reproduce the 
charge and current densities: 

p{r) = '^niip*{r)'^i{r) , (7) 

i 

i 

According to the variational principle and assuming a fixed set of the occupation numbers 
rii, we carry out the minimization of the Hohenberg-Kohn energy functional with respect to 
the single-electron orbitals '0j(r) subject to the constraint for the current density Q, i.e. the 
variation of auxiliary functional Q[$]. The minimization yields the following single-electron 
self-consistent equations: 

(^-iv + z;(r) + j dTj^^^+v,4p]iv) + ^^[d,,a{x)]^^Uv) = E,ij,{v) , (9) 



where the anti-commutator term, 

[d^, a{x)]. = d^{a{x)) + 2a(x)c?a; 



(10) 



is the additional imaginary potential arising directly from the constraint on electronic cur- 
rent (Appendix A). This equation is a particular case of the more general equation derived 
recently by Kosov and Greer ^. Relating a position gradient with the corresponding mo- 
mentum operator {p^ = —id^) we can rewrite the imaginary potential as the momentum- 
dependent operator: 

1 1 

(11) 



Since the operator |TT| is hermitian the Konh-Sham equation (jHl) is a hermitian eigenprob- 
lem. 

We have not specified yet the occupation numbers rij which enter the Kohn-Sham equa- 
tions (jni) via the density and the current density. To fix the set of occupation numbers, we 
need to know the number of electrons in the nanowire. Although we deal with the open 
quantum system where electrons can leave and enter the wire, we still can put a constraint 
on the total number of electrons because of time and space independence of the electronic 
current (i.e. steady state direct current). Then, the continuity equation tells us that the to- 
tal number of electrons does not depend upon time. If we assume that the charge neutrality 
is maintained after the establishment of the current currying state, the average number of 
electrons in the nanowire is the same as the number of electrons in the zero current wire. 
Therefore to obtain the occupation numbers we can minimize Q (jH)) with respect to rij sub- 
ject to the particle conserving constraint imposed with the additional Lagrange multiplier 
jj, (chemical potential): 



> . 



(12) 



The direct differentiation of the Hohen berg -Kohn energy is straitforward and results into 
the following expression (Janak theorem) [13]: 



drti 



vir 



dr- 



p(r) 



Vxc[p]{r) i>i{r) + 



drii 



v{v) + j c^r^^^ + v^,[p]{y)^ Mr) + c.c. 



The details of the differentiation of the constrained functional with respect to an occupation 
number rii are given in Appendix B. The differentiation leads to the following expression for 



the derivative: 



V'fc(r) + c.c. 



(13) 



Combining these two derivatives we arrive to the following expression for the derivative of 
the auxiliary functional: 

^ = 1 dr^;(r) (-iv + ^;(r) + J rfr^^ + t.,.,[p](r) + 1 [9., a(x)]^) ^.(r) + 
+C.C.] . (14) 



drii 



Given that ipi satisfy eq.Q it is easy to see that: 

- = eJ drij*{r)ij,{r) + ^ ^fc^fc^ / drijl{T)M^) = Ei . (15) 

From eq. (|T^ with the use of eq. (|TH|l we arrive to the inequality 

Y,{Ei - ^)5ni > (16) 

i 

which gives rise to the standard "build up principle" ( orbitals with Ei < fi have rij = 1 
and those with Ei > have Ui = 0) which is still hold within our approach: a maximum 
two electrons are put into current carrying orbitals in the order of increasing orbital energy. 
The approach can be readily extended to the finite temperature case replacing of the integer 
rij = or 1 occupation number by the Fermi-Dirac expression. The formal aspects of the 
DFT with direct current and non-zero temperature are discussed elsewhere. 

The Kohn-Sham equations with current are not yet in a form allowing for a solution 
to be found as the Lagrange multiplier a{x) is not known yet. The additional equation for 
the Lagrangian multiplier a{x) can be obtained if we require that the Kohn-Sham orbitals 
from the eq.ijH} yield the required current I{x): 

/" dydz{^P:{r)d,ij,{v)-^i{r)d,ij:{v)) = I{x) . (17) 

The system of nonlinear equations Q IT7|) are to be solved simultaneously and being com- 
bined together composes a set of Kohn-Sham equations for nanowire with direct current. 



IV. SELF-CONSISTENT ABSORBING-EMITTING BOUNDARY CONDITIONS 



Multiplying the eq. Q from the left on the njV'i(i')*, summing up over the Kohn-Sham 
orbitals and subtracting from the obtained equation its own complex conjugate yields the 
continuity equation: 

V-j(r) = 9,.(a(x)p(r)) . (18) 

After the integration of eq. (fT5j) over the {y,z} plane and substitution of the constraint for 
the current density eq.® into the equation the following first order differential equation for 
the Lagrange multiplier a{x) can be deduced: 

dx{a{x)py^{x) - I{x)) = , (19) 

where we have introduced the quantity 

Pyz{x) = / dydzp{r) . (20) 



The solution of the differential equation (fT^ depends upon the asymptotic behavior of 
the specified function for the current I{x). It is computationally convenient to assume the 
space-localized boundary conditions (/(x), Pyz{x) 0, when x ±oo). This space- localized 
boundary conditions yields zero integration constant for the differential equation (|T!I|l and 
we obtain the following expression for the Lagrange multiplier: 

With this choice of the Lagrange multiplier (PT|) we complete the Kohn-Sham equations with 
the current / in the self-consistent and closed form: 



I(x) d 



^,(r) = E.^,{v) . (22) 



_Pyz{xy dx_ 

Some aspects of eq.(|^ deserve special discussion. Suppose that the constrained current 
starts at the left system boundary L and completely absorbed at the right system boundary 
R: 

I{x) = 1 9{x- L)9{R - x) , (23) 

where / is the steady current through the system and 9 is the step function. Then one of 
the terms in the anticommutator imaginary potential becomes proportional to the derivative 
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of the step function and yields the singular imaginary 5-function potential on the system 
boundaries: 



This absorbing-emitting imaginary boundary potential is optimized due to its self- 
consistency. It depends upon the density and value of the input current and emits the 
current / at the left boundary and completely absorbs the same current at the right bound- 
ary. 

V. CONCLUSIONS 

In this paper, we have derived the set of Kohn-Sham equations for a nanowire with direct, 
steady current. Not restricted to the linear response, our approach uses a constrained 
minimization of the total energy with a subsidiary condition for the current density to 
formulate the Kohn-Sham equation with direct current. The subsidiary condition for the 
current density is maintained during the variation by a pointwise Lagrange multiplier. Being 
entirely based on the DFT, our method incudes electronic correlations into the transport 
calculations at the same level as they are taken into account into the corresponding exchange- 
correlation functional. 



li 



{S{x - L) - S{x - R)) . 
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APPENDIX A: FUNCTIONAL DERIVATIVES OF THE CONSTRAINT FUNC- 
TIONAL WITH RESPECT TO ORBITALS 



In this Appendix we compute the functional derivative of the constraint functional. The 
constraint functional has the following form 



with 



K[ijji,ijj*]= j dxa{x)(^j dydzj^ir)-!^ , (Al) 
^-(r) = 7: E (^*(r)a.^,(r) - iJ,{v)^,^P*{v)) (A2) 



2i 

The direct variation of the constraint functional with respect to the '0*(r) yields: 

k 



I \ S f 

2t 2t6ilj*{r)J ^ 

II 1 

= —a{x)d^i!i{v) + —d^{a{x)^pi{r)} = — [a{x),^^]^^p^{r) 



APPENDIX B: DIFFERENTIATION OF THE CONSTRAINED FUNCTIONAL 
WITH RESPECT TO AN OCCUPATION NUMBER 

In this appendix the technical details of the direct differentiation of the constraint func- 
tional ()A1|) is presented. 



d d f f 

Here we have introduced the new quantity which is the current through the orbital i inte- 
grated over the plane {y, z}: 
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Then 



drii 



A[ipi,ilj*] = J dxa{x)jix{x) + ^fc j dva{x) 

= / dxa{x)ji^{x) + nfc / dra(a;) 

— I dr 



5n,; 



9ni 



dm 



X o 

on, 

drii 



dxMr) - ^^d^^kir) ] + 



d^{a{x)7pk{r)) ^^d^{a{x)7pl{r)) ) = 



j dxa(x)jij;(x) + ^rik J dr ^^^J [g^, a(a;)]+-0fc(r) - c.c. 
^ y rfra(x) (^*(r)a^^i(r) - ipi{Y)d^i}l{Y)) + ^ j dr 



[d^,a{x)]j^ipk{Y) - c.c. 



1 

2i 



dvij*{v)[d^,a{x)]+il)i{r) + — 



^nfc / dr 



^|^[a,,a(x)]+V^,(r)-c.c. 



(B2) 
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